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ABSTRACT 

We present observational evidence of compressible magnetohydrodynamic wave modes propagating from 
the solar photosphere through to the base of the transition region in a solar magnetic pore. High cadence im¬ 
ages were obtained simultaneously across four wavelength bands using the Dunn Solar Telescope. Employing 
Eourier and wavelet techniques, sausage-mode oscillations displaying significant power were detected in both 
intensity and area fluctuations. The intensity and area fluctuations exhibit a range of periods from 181 — 412 s, 
with an average period ~290 s, consistent with the global p-mode spectrum. Intensity and area oscillations 
present in adjacent bandpasses were found to be out-of-phase with one another, displaying phase angles of 
6.12°, 5.82° and 15.97° between 4170 A continuum - G-band, G-band - Na I Di and Na I Di - Ca II K 
heights, respectively, reiterating the presence of upwardly-propagating sausage-mode waves. A phase relation¬ 
ship of ~0° between same-bandpass emission and area perturbations of the pore best categorises the waves 
as belonging to the ‘slow’ regime of a dispersion diagram. Theoretical calculations reveal that the waves are 
surface modes, with initial photospheric energies in excess of 35 000 W m~^. The wave energetics indicate a 
substantial decrease in energy with atmospheric height, confirming that magnetic pores are able to transport 
waves that exhibit appreciable energy damping, which may release considerable energy into the local chromo¬ 
spheric plasma. 

Subject headings: Magnetohydrodynamics (MHD) - Sun; Chromosphere - Sun; Oscillations - Sun; Photo¬ 
sphere 


1. INTRODUCTION 

An explanation of how the solar corona exhibits extra¬ 
ordinary temperatures exceeding 1 MK remains a key goal 


of solar research ffirdelyi||2004[ Klimchuk||2006[ Taroyan & 
Erdelyi|20091l. Of greater importance is how chromospheric 


plasma can maintain temperatures of up to 10 000 K, since 
heating this region demands an even g reater energy input 
to bal ance the extreme radiative losses ( jWithbroe & Noy^ 
\911). Atmospheric waves have long been proposed as the 
mechanism by which energy is transported from the solar sur¬ 
face to the upper atmosphere in order t o sustain its height¬ 
ened temperatures ( Schwarzschild|1948] l. The first verifiable 
evidence for oscillations i n the solar atmosphere was docu¬ 
mented by |Leighton] ( |1960] l, and in the years since this discov¬ 
ery, observations of global and localised wave activity have 
become more definitive as a result of much improved instru- 
mentatio n and processing techniques (see, e.g., the recent re¬ 
views by Banerjee et al.|2007| Mathioudakis et al.|2013[pess 
! these ol 


et al .|2015| l. Alongside these observational improvements, our 

erstandi 


understanding of solar oscillations has increased through the¬ 
oretical developments in the field of magnetohydrodynamics 


(MHD; e.g.. 

Roberts 2000[ 

Erdelyi|2008 

1 . From this, th 


e nature of the Sun’s atmosphere 


can be simplified as a region of magnetised plasma that can be 
manipulated through a number of physical mechanisms, in¬ 
cluding the ubiquitous convective flows present in the photo¬ 
sphere. MHD theory allows for an understanding of multiple 
wave modes in structured media, including sausage and kink 
modes, generated in the magnetically active regions of the so¬ 
lar atmosphere. The study of these waves in terms of their 
energetics and propagation through the solar plasma is vital 
in assessing the role oscillations play in atmospheric heating, 
given that the omnipresent solar acoustic (i.e., non-magnetic) 
waves have demonstrated insuffic ient energies to be the sole 
means of chromospheric h eating ( [Fossum & Carlsson]|2005[ 
[Bello Gonzalez et al.|2010l l. 

The various MHD wave modes that have been observed 
throughout the solar atmosphere can largely be categorised 
in terms of their compressibility. Weakly compressible MHD 
waves are distinguished by their inability to cause significant 
density perturbations within the plasma. Such modes iden¬ 
tified thus far are the kink wave, which has been identified 


both in the chromosphere (e.g., Kukhianic 

ze et al.|2006 Za- 

qarashvili et al.|2007jl and the corona (e.g.. 

Aschwanden et al. 

19991|Nakariakov et al.||1999UErdelyi & 1 

Laroyan|2008[ [Van 

Doorsselaere et al.||2008| to name but a few), alongside the 
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more e lusive and comple tely incompressible torsional Alfven 
wave (Jess et al. 2009|l. A consequence of the nearly in- 
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compressible nature of kink waves is their resistance to en¬ 
ergy dissipation unless lar ge Alfven speed gradients exist to 
help induce phase mixing ( Heyvaerts & Priest|1983|l, or pro¬ 
vide mod e conversion due to resonant coupling ( |Sakurai e~t| 
al.|1991[ l, making them the widespread focus of recent coro- 
nal heating studies. However, one drawback is that the pro¬ 
posed dissipation mechanisms for such waves exist on spatial 
scales below current observational resolution limits. As a re¬ 
sult, there has only been tentative observational signatures of 


fast kink mode wave damping in the chro mosphere ( |He et al.] 
2009 Kuridze et al.|2013t Morton|2014| l, though the mecha- 


nism remains largely undetermined and it is uncertain whether 
these waves contribute directly to chromospheric heating. 

Compressible wave modes are defined observationally by 
their ability to perturb the local plasma density. The greater 
capacity for mechanisms to extract energy from these waves, 
when compared to incompressible oscillations, continues to 
highlight their potential importance in providing heat energy 
to the outer solar atmosphere. The most observed compress¬ 
ible modes are the magneto-acoustic oscillations found in 
highly magnetic regions, such as sunspots, pores and mag¬ 
netic bright points. Inclined magnetic fields allow global solar 
acoustic waves to propagate further int o the atmosphere du e 
to a weakened net gravitational force (|Bel & Leroy [ 


Magneto-acoustic waves are thought to be associated witi 
other observed oscillations in the solar atmosphere. For ex¬ 
ample, running penumbral waves are observed to propagate 
outwards alon g the penumbral ma gnetic fields surrounding 
sunspots (e.g., Zirin & Stein||197^ . They are thought to be 
the chromospheric manifestation of photospheric magneto- 


acoustic waves as they propagate upwards (Bloomfield et al. 
2007[p^s et al.|2013|l, and have been shown to propagate up¬ 


wards tfoinTEecEromosphere with sufficient energy fl ux to 
contribute t o the localised heating of coronal plasma ( |Freij| 
et al. 2014| l. Compressible waves, in particular, are desir- 
able wave modes given that many forms of energy dissipation 
have been identifie d, including resona nt absorption (|Goossens 
& De Groof|2001 , turbulent mixin g ( [van Ballegooijen et al. 
201 l|l, viscosity ([BraginskiT] 1965 1 , rtiode con version iUlm- 
schneider et al. TTW and thermal conduction (|Ofman & As- 


chwanden|2002|l . 


Compressible MHD sausage modes are also of interest to 
the chromospheric/coronal heating debate. These waves are 
typically identified by periodic fluctuations in the intensity 
and/or area of magnetic structures suc h as pores, spicules or 
coronal loops ( |Edwin & Roberts|1983| l. Given the high spatial 
resolution required to observe such fractional changes in area, 
early studies of sausage mode waves were limited to larger- 


scale coronal loops (e.g., jAschwanden, Nakariakov & Mel 
nikov||2004)l, and it was only in recent years before the first 


lower atmospheric detection of sausag e mode waves in white 
light images of pores was uncovered ( |Dorotovic et al.|2008| l. 
As observing instrumentation has developed in capability, the 
ability to study sausage mode properti es at high spat i al and 
temporal resolutions has been evident. |Morton et al.| pOTT i 
provided the first in-depth study of sausage modes in photo- 
spheric pores. The identified area oscillations often did not 
have simultaneous intensity changes, indicating these waves 
did not possess large-amplitude wave power provided there 
was no twist in the magnetic field. H owever, these waves were 
ascertained to be fast mode waves ( |Moreels et al.|[20T3] l and 
their observed periods indicated that they may be driven by 
the globa l p-mode oscillations, i n contrast to the waves ob¬ 
served by|Dorotovic et al[](|2008[l, whose longer periods were 


postulated as evidence of the presence of a magneto-acoustic 
gr avity mode. Further pho tospheric analysis was conducted 
by Dorotovic et al. ( |2014 1 , where the investigation of sepa- 
rate pore and sunspot features identified the presence of both 
fast and slow sausage mode waves. These results showed that 
sausage modes form at the solar surface in a variety of dif¬ 
ferent configurations, and ha ve also been seen in fluctuat ions 
of spectro-polarimetric data ( |Fujimura & Tsuneta||2009| l, in¬ 
dicating a range of methods for photospheric detection and 
analysis. As the study of these waves at photospheric heights 
has developed, chromospheric detection s initially appeared to 
be more elusive. However, the work of [Morton et al.| ( |2012| l 
showed that sausage mode waves were in fact ubiquitous in 
fine-scale chromospheric magnetic structures, such as fibrils 
and mottles. These waves contained energy flux on the or¬ 
der of 11000 W m“^ and exhibited characteristics of leaky 
modes, suggesting that they were readily able to dissipate 
energy to the surrounding plasma. Given that the identified 
sausage-mode energy flux was greater than that of the simul¬ 
taneous kink modes also found in these structures, the work 
was important in validating their potential for energy trans- 
portation to the upper regions of the so lar atmosphere (e.g., 
[Jess et ar]|2012a[]Moreels et al.||2015a] l. However, what re¬ 
mains to be verified is a link between the photospheric and 
chromospheric observations of sausage modes, showing that 
the waves generated at the solar surface are indeed propagat¬ 
ing to higher atmospheric heights. 

Magnetic pores are the visible signature of concentrated 
magnetic fields that prohibit surface convection over a large 
volume. They are the smaller counterparts to sunspots, 
with diameters between 1- 6 Mm, and we aker magnetic field 
strengths up to ^ 1700 G ( |Sobotka|2003| l. There is an asso¬ 
ciated drop in the measured intensity, with the mean values 
of a pore being o n the order of 40% low er than the quies¬ 
cent solar surface ( jVerma & Denker|2014| l. Given their lack 
of penumbra, pores have often been consi dered (and mod- 
elled) as a simple single magnetic flux tube ([Simon & Weissj 
1970| Hurlburt & Rucklidge|2000 1 . However, observ ations of 


bright struc tures within pores, such as umbral dots (Sobotka 
jet al.||1999| ), could be an indication that pores may exhibit a 
m ore dynamic m agnetic configuration, such as that described 
by |Parker|(|1979|l. Their relatively small size has previously 
made detailed observations and analyses difficult, but with 
the advent of high resolution imaging instruments, they are 
now within scientific scope. This is advantageous since their 
small size leads to lessened internal forces when compared to 
sunspots, thus they are likely to be more dynamic and respond 
to bulk external forces more readily. 

Here, we present observations of upwardly propagating 
slow sausage mode waves in a solar magnetic pore using high 
resolution imaging and spectral imaging from the Dunn So¬ 
lar Telescope. We also provide the first intensive analyses of 
sausage mode energies with atmospheric height, from their 
formation in the photosphere through to high chromospheric 
locations. 

2. OBSERVATIONS AND DATA PROCESSING 

The data presented here is an observational sequence ob¬ 
tained during 19:27 — 20:02 UT on 2013 March 6 with the 
Dunn Solar Telescope (DST) at Sacramento Peak, New Mex¬ 
ico. The telescope was focussed on the active region NOAA 
11683, positioned at heliocentric co-ordinates (113", —170"), 
or S17.3W07.0 in the conventional heliographic co-ordinate 
system, with seeing conditions remaining very good through- 
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Fig. 1 .— The left column consists of co-spatial images of the full ROSA/CSUNcam field-of-view, stacked from the photosphere through to the chromosphere. 
From bottom to top, the images represent the vertical magnetic field strength, Bz, determined from HMI vector magnetograms with the scale saturated at 
±1000 G to aid clarity, magnetic field inclination angles, where 0° and 180° represent fields outwardly and inwardly normal to the solar surface, ROSA 4170 A 
continuum, G-band, Na I Di and Ca II K images. A white dashed line interconnects the pore between bandpasses, while the solid boxes define the sub-fields 
displayed in the right-hand column. Here, the pore is shown in detail for each imaging bandpass, with the binary map pixels shown at each height using a red 
contour. 


out the entire 35 minute observing rvin. The Rapid Oscil- 
lations in the Solar Atmosphere (ROSA; |Jess et al.||2010i 
camera system was used to image a 115" x 115" portion of 
the solar disk through continuum (4170 A; 52 A FWHM), G- 
band (4305.5 A; 9.2 A FWHM) and Na I Di core (5895.9 A; 
0.17 A FWHM) filters, each at a spatial sampling of 0712 
per pixel. A Ca II K core (3933.7 A; 1.0 A FWHM) filter 
was also used in conjunction with a new camera addition 
to the DST’s imaging suite. The California State University 
Northridge camera (CSUNcam) is an iXon X3 DU-897-EX[|] 
model manufactured by Andor Technology, boasting a back- 
illuminated, 512x512 pixeP electron-multiplying CCD. The 
detector, triggering and readout architectures are identical to 

* Full specifications available at http://www.andor.com 


that of the exist ing Hydrogen-alp ha Rapid Dynamics cam¬ 
era (HARDcam; [less et al.|2012c), but incorporates an addi¬ 


tional window coating designed to improve sensitivity in the 
blue portion of the optical spectrum, helping quantum effi¬ 
ciencies exceed 70% at 4000 A. To ensure the solar coverage 
of CSUNcam was similar to that of ROSA, a plate scale of 
07230 per pixel was chosen to provide a 118" x 118" field- 
of-view. 

In addition to ROSA and CSUNcam observations, the In¬ 


terferometric Bldimensional Spectrometer (IBIS; Cavallini 
2006|l was used to simultaneously sample the Ca 11 absorp 


tion profile at 8542.12 A with a spatial sampling of 07097 per 
pixel. IBIS employed 11 discreet wavelength steps, each with 
10 exposures per wavelength to assist with image reconstruc¬ 
tion, providing a complete scan cadence of 17.673 s. A white- 
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light camera, synchronised with the IBIS feed, was used to 
enable the processing and destretching of all narrowband im¬ 
ages. 

To improve the clari ty of the images , techniques of high- 
order adaptiye_optic^ Rimmek 2QQ4\ , speckle rec onstruc- 
^ " ‘Jess et 


tion (|W6ger et al.||2008[ ) and Fourier co-alignment (jji 
al.|2007|l were employed. Utilising 64 —1 speckle re 


restora¬ 


tions, the resulting cadences for continuum, G-band, Na I Di 
and Ca II K image sequences were 2.11s, 2.11s, 4.22 s and 
2.33 s, respectively. The IBIS images were processed with a 
10 —> 1 restoration, and a blueshift correction was applied to 
ac count for the use of classical etalon mountings (|Cauzzi el 


ah]|2008 1. The Na I Di and Ca II K image sequences were 


then interpolated on to a constant 2.11s time grid to allow 
for easy comparisons between bandpasses, with the Ca II K 
images also spatially resampled to replicate the plate scale of 
the ROSA cameras. The end result is spatially and tempo¬ 
rally coaligned image sequences spanning the full 35 minute 
duration of the data set. 


The Helioseismic and Magnetic Imager (HMI; Schou et ah 
2012|l onboard the Solar Dynamics Observatory (SDO; Pes- 


'nell et all2012|l was utilised to provide simultaneous vector 


magnetograms of active region NOAA 11683, in addition to 
contextual HMI continuum images for the purposes of co- 
aligning to the ground-based data sets. The Milne-Eddington 
vector magnetograms were obtained with a two-pixel spatial 
resolution of I'.'O and a cadence of 720 s. The HMI data were 
processed using the standard hini_prep IDL routine, which 
includes the removal of energetic particle hits. Next, 200" 
X 200" sub-helds were extracted from the processed images, 
with a central pointing approximately equal to that of the 
ground-based data. Using the HMI continuum image to de¬ 
fine absolute solar co-ordinates, the ROSA and CSUNcam ob¬ 
servations were subjected to Fourier-based cross-correlation 
techniques to provide sub-pixel inter-bandpass co-alignment. 
To do this, the plate scales of the ROSA and CSUNcam ob¬ 
servations were hrst degraded to match that of the HMI con¬ 
tinuum imag^ Then, calculations of squared mean abso¬ 
lute deviations were performed between the datasets, with the 
ground-based images shifted to best align with the HMI ref¬ 
erence image. Following the implementation of co-alignment 
techniques, the maximum x- and y-displacements were both 
less than one tenth of an HMI pixel, or 0"05 («36 km). Sam¬ 
ple ROSA, CSUNcam and HMI images are displayed in Fig¬ 
ure [U 


3. ANALYSIS AND DISCUSSION 

Within our held-of-view there was a pronounced pore struc¬ 
ture, as highlighted by the solid white boxes in Figure that 
lasted for the full duration of the data set. The pore was a 
well-dehned dark discontinuity from the surrounding plasma. 
Thus, to isolate the pore for further study, we employed in¬ 
tensity thresholding techniques allowing the perimeter of the 
pore to be dehned in relation to a decrease below the mean 
quiescent intensity. A quiet region, free from magnetic struc¬ 
tures such as plage and bright points, was used to calculate the 
background mean intensity, /ave, and the resulting standard 
deviation, a. Next, the pore was dehned in each bandpass as 
pixels displaying intensity values less than (lave —where 
N is equal to 2.2, 2.2, 2.2 and 2.1 for the 4170 A continuum, 
G-band, Na I Di and Ca II K bandpasses, respectively. While 

^ Data analysis was performed on full-resolution (i.e. non-degraded) im¬ 
age sequences. 


the value of N may differ between adjacent bandpasses as a 
result of variable height-localised contrast ratios, once dehned 
it remained hxed for each bandpass over the entire duration of 
the time series to ensure consistency when measuring time- 
dependent variables. 

By employing the intensity threshold values, pore time se¬ 
ries for each bandpass were generated by isolating those pix¬ 
els which lay at or below the specihc threshold value. This 
created a 988 (35 minute duration with a 2.11 s cadence) im¬ 
age sequence for each bandpass with all solar structures ex¬ 
cept for the pore masked out. Then, the pore areas were time- 
averaged in each bandpass to establish the mean pore cover¬ 
age at individual atmospheric heights. The time-averaged area 
of the pore dehned by the intensity thresholds corresponded 
to 7.02 ±0.11 MmA 7.18 ± 0.11 Mm^, 7.99 ± 0.13Mm2 and 
8.12 ± 0.19Mm^ for the 4170 A continuum, G-band, Na I Di 
and Ca II K bandpasses, respectively, indicating that the pore 
presented here expands minimally with height (^15% area 
change between continuum and Ca II K layers). This expan¬ 
sion with height is expected due to the decrease in plasma 
pressure, although the fr actional change in are a is smaller than 
typically expected (e.g., [Solanki et al.|1999] ). This may be a 
direct consequence of the strong magnetic helds present in the 
region of the pore at photospheric heights (left column of Fig¬ 
ure [T]), which prevent rapid expansion as it extends upwards 
from the continuum level, or an unresolved twist in the held 
lines that would result in reduced expansion. With the pore 
perimeter well dehned in each bandpass, it became possible 
to investigate and compare propagating and height-localised 
Huctuations in intensity and area, as will be discussed in the 
following subsections. 


3.1. Oscillations in the Pore Intensity 


To investigate whether the pore structure contained stand¬ 
ing and/or propagating magneto-acoustic wave trains we stud¬ 
ied co-spatial intensity Huctuations between adjacent band¬ 


passes using the wavelet analysis techniques described by Jess 


|et al.| ( |2012b| l, and hrst introduced by Torrence & Compo 
( |1998[ )! First, however, it was important to isolate pixels that 
remained part of the pore structure at all points in time and 
height. Only with these pixels determined can wavetrains that 
propagate through all available bandpasses be identified and 


statistically studied. Given that the 4170 A continuum chan¬ 
nel represented the smallest spatial coverage of the pore, it 
was used to determine which pixels are continually contained 
within the pore boundaries. The binary map was formed by 


co-adding all 4170 A continuum pore perimeter maps, where 
values of ‘ 1 ’ indicated part of the pore structure and values of 
‘0’ represented external non-pore areas, normalising by the 
image count (988), then setting all pixels with normalised 
intensities <1 to zero. Thus, only pixels that are contained 
within the pore perimeter for all steps in time are assigned 
a value of ‘1’ in the hnal binary map, as indicated by the 
red contours in the right-hand panels of Figure [T] The pix¬ 
els contained within the 4170 A continuum binary map were 
then compared to the time series for other bandpasses to verify 
their omnipresence throughout time and atmospheric height. 
The completed binary map contained 206 non-zero pixels, 
and each pore image sequence was subsequently multiplied 
by the binary map to form a data set that could be compared 
directly between all atmospheric heights and time steps. 


Intensity time series for the 4170 A continuum and G-band 
data sets were extracted for those pixels which lay within 
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Power Spectrum of Calcium Pore Intensities 




Fig. 2. — The top panel is the sum of individual Fourier power spectra derived within the pore structure at Ca IIK wavelengths. The Fourier power is displayed in 
arbitrary units, while the oscillatory periods are plotted on a logarithmic scale to better display oscillations synonymous with the longer-period p-mode spectrum. 
The lower panels are histograms of the phase angle values from the pixel-by-pixel intensity analyses for the 4170 A continuum - G-band, G-band - Na I Di and 
Na I Di - Ca II K bandpasses. For each inter-bandpass gap, a peak can be defined representing the most prominent phase angle observed. 


the confines of the binary map. The time series were then 
detrended by a first-order polynomial to remove long term 
variations in intensity and/or light levels, and normalised to 
their subsequent mean. Wavelet analysis techniques were ap¬ 
plied to the time series to identify the presence of periodic 
signals. Strong oscillatory power was detected in both band- 
passes, with periodicities predominantly in the range of 3 to 
5 minutes, which overlaps well with the typical solar p-mode 


Spectrum (|Lites et al.|1982 

1 . Next, the phase difference anal- 

ysis techniques detailed by 

Jess et al.|(|2012b|l were employed 


to investigate whether the oscillations were detectable co- 
spatially in the neighbouring bandpasses, and whether there 
was an inherent phase shift. Here, a positive phase angle in¬ 
dicates that a wave is first observed at a lower atmospheric 
height (i.e., the wave is propagating upwards from the 4170 A 
continuum to the G-band), with criteria defined to ensure con¬ 
fidence in the final results. Firstly, only waves that exhibited 
a cross-correlation coefficient greater than 50% between the 
bandpasses were considered in an attempt to ensure the os¬ 
cillatory signature corresponded to the same wave. Secondly, 
only waves which displayed a lifetime greater than ^/2P in 
each bandpass were included, where P is the period of the 
oscillation in seconds, to ensure these signatures are periodic 
rather than lone spikes in the time series. These rigorous cri¬ 
teria were implemented on each extracted time series, with 
approximately 1700 suitable periodicities identified in the pe¬ 
riod range 210 — 412 s, with a mean value of 290 ± 31 s. The 
derived individual phase angles between these detected os¬ 
cillations were collected into a histogram employing 1° bins 


(as detailed in Figure |^. A best-fit Gaussian was applied, 
with the centroided value identifying the most commonly oc¬ 
curring phase angle. In the case of 4170 A continuum - G- 
band, the dominant angle found was «6.12 ± 4.6°, indicat¬ 
ing upward propagation between the formation heights of the 
4170 A continuum and G-band filtergrams. Subsequently, the 
same methodology was applied to the G-band - Na I Di and 
Na I Di - Ca II K time series. In the case of G-band-Na I Di, 
periodicities in the range of 210 — 412 s were identified (sim¬ 
ilar to that of the 4170 A continuum - G-band), with a promi¬ 
nent phase angle of «5.82±3.62°. The Na I Di - Ca II K out¬ 
put displayed a small decrease in the range of detected peri¬ 
odicities, with periods in the range of 181 — 362 s being found 
to display a dominant phase angle wlS.Q? ± 8.73°. The up¬ 
per panel of Figure [^displays a Fourier transform of Ca II K 
intensities synonymous with all pixels contained within the 
pore structure, here demonstrating significant power at peri¬ 
odicities in the range of 200 - 300 s. These relationships indi¬ 
cate the prevalence of upwardly propagating waves within the 
pore. 

From inspection of Figure S the majority of the derived 
phase angles are positive, witli75% of pixels showing a pos¬ 


itive phase between 4170 A continuum - G-band, 70% be¬ 
tween G-band - Na I Di, and 78% between Na I Di - Ca II K. 
The percentage of upwardly propagati ng waves de t ermined 
here compares well with the work of |Jess et al.| (|2012b||, 
which ascertained that 73% of waves propagating between 
the 4170 A continuum and G-band channels within kiloGauss 
magnetic elements were directed upwards. The existence of 
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downwardly propagating waves is to be expected, and ulti¬ 
mately may be a consequence of reflections occurring at the 
chromospheric and transition region boundaries. In addition, 
the formation heights of the Na I Di and Ca II K flltergrams 
are difficult to quantify, with contribution functions being 
comprised of emission generated over a range of heights and 
opacities. As such, the phase angles calculated between these 
bands will have an intrinsic degree of scatter as a result of non- 
fixed formations heights and varying opacities in the vicinity 
of the highly magnetic pore. However, even with such effects, 
each inter-bandpass jump still demonstrates in excess of 70% 
upwardly propagating wave signatures. 

The phase speed of the waves can be estimated by combin¬ 
ing the observed phase angle alongside the formation heights 


of the 4170 A cont inuum and G-band bandpasses derived by 
Jess et ari(2012b]l. While the average separation between 


formation heights was found to be ^75 km, structures re¬ 
siding within more magnetic environments (such as inter¬ 
granular lanes and magnetic bright points), were found to 
display larger separation heights on the order of ~100 km. 
While more difficult to constrain, the formation heights of 
the Na I Di and Ca II K flltergr ams can be est imated as 
~400 kr n (jSimoniello et al.||2008l l and ^800 km ( |Beebe &| 
|Johnson||l969| l, respectively. Employing our periodicity and 
phase delay measurements from the 4170 A continuum to G- 


band, along with the traversed distance extracted from Jess et 


al. (|2012b|l, enables us to calculate the resulting phase veloc¬ 


ities of the observed waves. Assuming the prominent phase 
angle of 6.12°, found between 4170 A continuum and G-band 
flltergrams, represents the phase velocity of the waves with a 
mean period of 290 s, the wave travel time can then be cal¬ 
culated as 6.12°/360° x 290 s = 4.92 s. Thus, an average 
propagation distance of ~75 km leads to an estimated wave 
speed of ~15 km s“^. At first glance, this value is in excess of 
what is expe cted for an isolated wave propagating in the solar 
photosphere (jEdwin & Roberts|l983|l. However, investigation 
of the line-of-sight Doppler and bi-sector velocities obtained 
with IBIS revealed a large bulk upflow within the confines 
of the pore. Eigure [^displays the spatially averaged Doppler 
and bi-sector velocities within the pore at various positions in 
the line profile. As can be seen from the bi-sector velocities, 
which best represent lower atmospheric heights, a global up¬ 
flow with an average velocity of ^8 kms“^ is present, with 
a peak value approaching 14 kms“^ recorded. The observed 
decrease in the plasma upflow velocities towards the Ca II line 


core re mains consistent with the previous work of Cho et al. 
(|2013|). With the addition of a bulk plasma upflow, the initial 
speed calculated can be considered as a Doppler-shifted phase 
speed, which incorporates both the phase velocity of the wave 
and the bulk motion of the transfer medium. Therefore, an 
estimate of the true phase speed of the waves can then be de¬ 
rived by subtracting the line-of-sight measurements from the 
derived velocity, thus producing a phase speed on the order 
of 5 kms“^ at lower atmospheric heights, which is similar 
to, if not below, the adiabatic sound spe ed predicted for these 
locations (see, e.g., Jess et ari2012b| l. Thus, the effect of 
the waves being superimposed on top of a simultaneous bulk 
upflow, alongside the p otential effects of a twis ted magnetic 
held configuration (e.g.,|Erdelyi & Eedun|2006|l may account 
for the large measured Eioppler-shifted phase velocity. 


3.2. Oscillations in the Pore Area 


Observed LOS Doppler Velocities within the Pore 



0 500 1000 1500 2000 

Time (s) 


Fig. 3.— The average line-of-sight velocities within the magnetic pore, 
measured from the IBIS 8542 A Ca II data. The blue line represents the ve¬ 
locity shifts in the line core, with the red and green lines denoting the velocity 
shifts at 100 mA and 150 mA, respectively, representing lower heights in the 
atmosphere away from the chromospheric line core. The line core average 
velocity is 3.27 kms“^, with the 100mA and 150mA bi-sectors exhibiting 
average velocities equal to 5.93 km s“^ and 8.01 kms“^, respectively. 


Upwardly propagating intensity oscillations ma y be a direct 
signature of magneto-acoustic wave modes (e.g.,|De Pontieu 
et al.|2004|Bess et al.|2012^ Ereij et al.|2014 1 . ITierefore, in 
order to verify whether the detected oscillation 


nations also have the 


observational signatures of sausage mode waves, it was nec¬ 
essary to consider whether simultaneous periodic changes in 
the area of the pore existed. When compared to the pixel-by- 
pixel analysis employed for the intensity oscillations, which 
provided 206 individual pixel measurements per time step, 
the examination of periodic area changes has statistical limi¬ 
tations due to only one area measurement being possible per 
instant in time (rather than 206 individual intensity measure¬ 
ments). Nevertheless, the extraction of phase delays between 
pore area sizes as a function of atmospheric heigh t was per¬ 
formed identically to that described in Section [3T] 

The analyses identified area oscillations of similar period 
ranges (181 — 412 s) to that of the intensity oscillations, thus 
indicating the presence of compressional waves that are able 
to modify the cross-sectional area with time. The results of the 
phase analysis on neighbouring bandpasses again provided a 
consistently positive phase angle, reiterating the presence of 
upwardly propagating compressive waves. Specifically, the 


phase delays found in area oscillations between the 4170 A 
continuum - G-band, G-band - Na I Di and Na I Di - Ca IIK 
bandpasses were 6.35±2.16°, 5.72±6.7° and 16.76±13.3°, 
respectively. These results are comparable to the peak values 
found for the pixel-by-pixel intensity oscillations (see Eig- 
ure|^. Thus, the identification of synchronous intensity and 
area oscillations, alongside the remarkable similarity between 
the phase delays for both sets of perturbations, provides evi¬ 
dence for the presence of upwardly propagating sausage mode 
waves within the pore. 


3.3. Determining the Wave Mode 

Employing the phase analysis techniques independently to 
both intensity and area fluctuations as a function of atmo¬ 
spheric height has confirmed the presence of upwardly prop¬ 
agating sausage mode oscillations. However, these oscilla- 
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tions may fall within the ‘slow’ or ‘fast’ sections of the dis¬ 
persion relation. The categorisation of waves in this manner 
allows limita tions to be imposed upon properties such as the 
phase speed (Edwin & Roberts|1983]|Evans & Roberts|1990[ 


Erdelyi & Fedun||2010|l and dissipative potential through 

mechanisms such as thermal conduction ( |Van Doorsselaere] 
et al. 201 la[ ). Therefore, in order to categorise whether 


the observed sausage mode oscillations are ‘fast’ or ‘slow’, 
the phase relationships between localised intensity and La- 
grangian area oscillations within each respective bandpass 
was investigated. Theory predicts an in-phase relationship be¬ 
tween intensity and area for slow-mode waves, while a phase 
lag of ^180° sugges ts the presence of fast-mode oscillations 
( |Moreels et al.|20f3] l. 


The criteria used in sections 3.1 and 3.2 were implemented 
to calculate the relevant phase relationships for the four band 
passes. Prominent phase delays between intensity and area of 
17 ± 21°, 19 ± 16°, 11 ± 3° and 21 ± 20° were found for the 
4170 A continuum, G-band, Na I Di and Ca II K bandpasses, 
respectively, with the area consistently leading the intensity 


fluc tuatio ns in each bandpass. In accordance with Moreels et 
[aE] ( |2013| ), this suggests a dominant presence of slow mode 
waves. 

Sausage mode waves have already been shown to exhibit 
strong velocity components along the tube axis (|Fujimura 
& Tsuneta||2009|l. As such, oscillatory patterns in line-of- 
sight (LOS) Doppler velocities can be employed as a diag¬ 
nostic tool when comparing with simultaneous intensity fluc¬ 
tuations. Phase relationships between the LOS Doppler ve- 


locity and the Eulerian intensity are identified in Moreels & 
Van Doorsselaere (2013|l. It is shown that slow propagat- 


ing modes will demonstrate in-phase behaviour. The LOS 
velocity can be determined from the Ca II data from IBIS, 
where a series of velocity bi-sectors were calculated. Such bi¬ 
sector measurements sample the plasma velocities at locations 
in the profile wings (i.e., away from the chromospheric line 
core), and are therefore analogous to velocity measurements 
at lower heights in the solar atmosphere. Thus, the LOS val¬ 
ues calculated from the wings of the Ca II IBIS scans will be 


formed higher than the pure solar continuum (i.e., the 4170 A 
ROSA continuum), yet lower than the chromospheric forma¬ 
tion heights of the Ca II line core, including those found in 
the Ca II K ROSA filtergrams. Therefore, in the case of the 
present data, the velocity bi-sectors are most compatible with 
the Na I Di observational time series. Oscillations found in 
the 150 mA bi-sector velocities are observed to have a phase 
relationship of «0.4 ± 0.1° with the Eulerian intensity oscil¬ 
lations derived from the Na I Di filtergrams. 

The approximately in-phase relationship between the LOS 
Doppler velocities and the Eulerian intensities provides confi¬ 
dence in a slow mode interpretation. The slight discrepancies 
found between observation and theory in terms of the area and 
intensity phase angles (i.e., 17°, 19°, 11° and 21°; see above) 
may be an indication that there is a superposition of both fast 
and slow sausage modes within the pore (albeit dominated by 
the slow modes). It may also be the result of varying opac¬ 
ity across the diameter of the pore, hence inducing variations 
in the pixel-by-pixel formation heights sampled by the filter- 
grams. This phase la g could also be caused by diss ipative 
effects, as observed by |Van Doorsselaere et al.| ( f201 lb| l. How¬ 
ever, due to the clear in-phase relationship between the LOS 
Doppler velocities and the Eulerian intensities, we are able 
to best interpret the observed waves as propagating slow ax- 


isymmetric modes. 


3.4. Assessing the Propagating Energy Flux of the Waves 

Using a newly developed fram ework for evaluating t he en¬ 
ergy flux of axisymmetric waves (|Moreels et al.|2015b|l, a de¬ 
tailed analysis of the energy of the observed MHD waves at 
each atmospheric height can be deduced for the first time. Be¬ 
fore the energy analyses can be undertaken, features of both 
the pore, its surroundings and individual sausage oscillations 
must be estimated as inputs to the energy calculations. 


3.4.1. Deriving Plasma Parameters 

In order to investigate the energy content of the waves, lo¬ 
cal plasma parameters must be identified that can then be inte¬ 
grated into a similar model atmosphere. Empl oying the Very 
Fast I nversion of the Stokes Vector (VFISV; |Borrero et al.| 
201 l| l algorithm, the line-of-sight magnetograms provided by 
the HMI instrument were decomposed to define a magnetic 
component perpendicular the solar surface (i.e., Bf). The 
formation height of the HMI ma gnetograms corresponds to 
~300 km above the photosphere ([Norton et al.||2006l |Fleck| 
et al.||2011|), with the pore revealing itself as a positive po¬ 
larity wim an area-averaged vertical magnetic field strength 
wbOO G, peaking at nearly 1500 G in the centra l region of the 
pore (l ower-left panel in Figure [^. Next, the |Maltby et al.| 
|T986|) ‘hot’ umbral model ‘L’ was scaled so the magnetic 
lield strength of the model (at a height of ^300 km) matched 
that of our HMI observations. The ‘hot’ umbral model was 
chosen as the best approximation to the pore, given that it is 
significantly smaller than a typical sunspot structure, and as 
a result is likely to be marginally hotter at its core due to the 
fact that convective proce sses will not be as severel y inhib- 
ited (jSobotka et al.|1999|). Utilising the (non-scaled) Maltby| 
et al.|(|1986|l density values, alongside the scaled magnetic- 
held strengths, enabled the calculation of the internal Alfven 
speed, Ua, as a function of atmospheric height using the rela¬ 
tion, 

Bz 

Va = —^ , (1) 

where B is the magnetic field strength, /tq is the magnetic 
permeability and p is the density of the local plasma. The 
solid black lines in Figure |^displays the plasma densities, ki¬ 
netic temperatures, scaled magnetic field strengths and Alfven 
speeds as a function of atmospheric height. It is clear that for 
lower atmospheric heights up to ^700 km, the Alfven speed 
determined is ~ 12 km s “ ^. 

It was also necessary to employ model atmosphere pa¬ 
rameters, corresponding to solar plasma conditions outside 
the magn etic pore, that could be directly compared with 
the scaled [Maltby et al.| (jl98^ values. Examination of the 
bottom-left panel of Figure |T| re veals that the pore exterior is 
still magnetic in nature, witn field str engths of the order of a 
few hundred Gauss. As a result, the Vernazza et al.| ( |1981| l 
model ‘D’ atmosphere was chosen to best represent the exte¬ 
rior of the pore structure. The VAL-D parameters are consis¬ 
tent with magnetic network elements, and thus best represent 
the nature of the solar plasma immediately outside the pore 
structure. However, as the atmosp heri c height scales differ 
between the [Vernazza et al.| ([1981 1 and [Maltby et al.[ ([1986 1 
models, each set of parameter measurements were interpo¬ 
lated on to a constant height separation grid equal to 100 km. 
The resulting VAL-D plasma densities, kinetic temperatures. 
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Grant et al. 


Plasma Density Kinetic Temperature 
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Alfven Speed 



Fig. 4.— The plasma density (upper-left), kinetic temperature (upper-right), vertical magnetic field strength (Szi lower-left) and deriv ed Alfven speed (lo wer- 
right), estimated as a function of height above the solar photosphere for internal (solid line; deduced by scaling the ‘hot’ umbral model of|Maltby et al.|1986f and 
external (dashed line; estimated directly from the VAL-D model of|Vernazza et al.|198H pore configurations. 


magnetic field strengths and derived Alfven speeds are dis¬ 
played in Figure |4| using a dashed lin e. Of notable compari¬ 
son to the scaled [Maltby et aL] ! 1986| l values are the tempera¬ 
tures, magnetic field strengths and Alfven speeds, all of which 
are smaller within the lower solar atmosphere, thus remaining 
consistent with our ROSA/CSUNcam imaging and HMI mag¬ 
netogram observations. 

Although these models provide one of the best representa¬ 
tions of the pore and its surroundings currently available, their 
accuracy of the true conditions observed may be limited. In 
particular, the modelled magnetic field is scaled through ob¬ 
servations at a single isolated height, leading to potential un¬ 
certainties as a result of complex magnetic structuring. This 
may produce an over-estimation of magnetic held strength in 
the lower atmosphere, especially considering our derived held 
strength at the solar surface is on the order of 3000 G. Further¬ 
more, a typical sunspot model will incorporate atmospheric 
expansion greater than the 15% we observed for the pore. As 
a result, this may lead to the modelled pore exhibiting greater 
plasma densities at higher atmospheric heights compared to 
real life, and as such, the internal Alfven speed calculated here 
may be large in relation to typical chromospheric channels. 
The lack of measured expansion may also lead to larger longi¬ 
tudinal magnetic held strengths at chromospheric heights than 
predicted by the model atmosphere as a result of magnetic Hux 
conservation. However, given that the pore structure is not ob¬ 
servable at higher heights in simultaneous SDO imaging, this 
suggests that the magnetic held may have fanned out into the 
surrounding magnetic canopy, as one would expect for larger 
sunspot structures at chromospheric and coronal heights. As a 


consequence, such fanning will lead to a progressive decrease 
in the longitudinal magnetic held strength with atmospheric 
height, in agreement with the Maltby model (see the lower- 
left panel of Figure |^. In future, additional observational 
steps can be taken to mitigate the issues highlighed here. In 
particular, complimentary spectro-polarimetric observations, 
using absorption lines such as Fe I 6302 A, would allow for 
an accurate measurement of the magnetic held at a different 
atmospheric height, hence better constraining the magnitudes 
of the magnetic helds. While we openly discuss any potential 
limitations of the model parameters derived here, they still re¬ 
main one of the best approximations currently available. 

3.4.2. Identifying and Characterising Individual Perturbations 

Thus far, the axisymmetric oscillations have been con¬ 
sidered as a superposition of multiple individual periodici- 
ties. However, in ord er to implement the energy analysis of 
[Moreels et al.| ( |2015b| l, the oscillatory properties of individual 
periodicities must be extracted from the time series. It was 
therefore necessary to isolate linear oscillatory perturbations 
in order for the energy analysis to be valid, but also to isolate 
waves that could be identihed over a narrow period range in 
order to limit the associated uncertainties in the frequency do¬ 
main. Periodicities were investigated by hltering the intensity 
and area time series in the frequency domain through use of 
Fourier transforms. Two periods in particular were identified, 
a 210 s sausage mode that was present for the first 844 s of 
the observations, and a 290 s wave that was present for the 
entire 35 min duration of the time series (see, e.g., the upper 
panel of Figure |^. Both wave modes exhibited clear sinu- 
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TABLE 1 

Properties of the upwardly propagating 210s and 290s oscillations 


Parameters 

Continuum 

G-Band 

NalDi 

Ca II K 


210s 

290s 

210s 

290s 

210s 

290s 

210s 

290s 

Intensity Perturbation (%) 

3.6 ± 2.5 

3.4 ± 2.4 

3.6 ± 2.5 

2.8 ± 1.9 

2.95 ± 2.09 

2.2 ± 1.5 

5.53 ± 3.91 

5.6 ±4 

Area Perturbation (%) 

1.53 ± 1.08 

2.1± 1.48 

2.92 ± 2.07 

2.1 ± 1.48 

3.25 ± 2.29 

3.25 ± 1.62 

4.75 ± 3.36 

5.8 ±4.1 

Phase Speed (km/s) 

4.03 ± 0.66 

2.5 ± 1 

2.26 ± 0.83 

1.9 ±0.68 

2.09 ± 0.89 

2.08 ± 0.9 

3 ± 1.5 

2.68 ± 1.3 

Energy Flux (kW/m^) 

54.5 ± 26.3 

38.4 ±21.4 

70.1 ± 35.8 

34.1 ± 13.5 

1.99 ±0.8 

2.48 ± 1.03 

0.12 ±0.08 

0.13 ± 0.08 


soidal behaviour, allowing area and intensity amplitudes to be 
easily extracted as a percentage fluctuation above the mean, 
as detailed in Table 1. 

An important input to the energy calculations is the phase 
speed of the waves as they p ropagate upwards, and the meth¬ 
ods of|Moreels et al.|(|2015all were used to calculate the phase 
speed at each atmospheric height. The methodology utilises 
a uniform flux tube as an equilibrium model, and combined 
with linear MHD theory, establishes a link between the phase 
speed of an oscillation and the amplitude of the total intensity 
and area fluctuations. A dimensionless amplitude ratio, A is 
defined as, 

where SI is the amplitude of the intensity perturbation, Ig is 
the mean intensity, SS is the amplitude of the area perturba¬ 
tion and So is the mean area. The phase speed is therefore 
defined as. 



±A-1 


±A — 1 -f (7 — l){hv/kBT) 


(3) 


where uj/k is the phase speed, 7 is the ratio of specific heats, 
h is the Planck constant, ly is the frequency at which the ob¬ 
servation was taken, fcs is the Boltzmann constant, and T is 
temperature within the pore. The sound speed, eg, is defined 


cs = 


hRT 


(4) 


where 7 is the ratio of specific heats, R is the gas constant, T 
is the temperature and /r is the mean molecular weight. How¬ 
ever, equation is, strictly speaking, only valid for contin¬ 
uum emission and not for specific emission/absorption lines. 
Therefore equation should only be applied to the continua 
found in the 4170 A^ G-band, Na I Di and Ca II K filters. 
While significant numbers of absorption features are present 
in both the 4170 A and G-band filter response curves, the large 
spectral coverage of the filters means that they will predom¬ 
inantly be comprised of solar continua. For the Na I Di and 
Ca II K image sequences we need a different inversion for¬ 
mula. Let us assume that for any optically thick absorption 
line we can write the intensity in terms of the local density 
and temperature, i.e., I = F{p,T) where F" is a smooth 
function. Following the same methodology as IMoreels et al. 
( |20T3i|2015al l we find. 


UJ 




±A- 1 


±A-1 


£0 I / 1 N ^ 

Fo dp ^ y I Fo dT 


(5) 


To respectively. The indices ‘zero’ stand for the equilibrium 
parameters, i.e., Fq = F{pq,To). Indeed, when we use 
F{p,T) = (2/ii/^/c^) exp(—/ii//fcBT) we recover equation 
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Unfortunately, the use of equation requires a detailed 
knowledge of the precise formation heights of the Na I Di 
and Ca II K spectral lines. Examining other lower atmo- 
spheric spectral features, incl uding the Ha and Mg II lines, 
Leenaarts et al.| (|2012| |20T3]l demonstrated the difficulty in 
accurately linking the observed intensity to the localised den¬ 
sity/temperature. The authors claim that due to the chromo¬ 
sphere not being in local thermodynamic equilibrium, it is 
therefore difficult to isolate and qu antify the local thermal 
properties of the plasma. However, |Leenaarts et al.| (|2013 1 
revealed that the temperature at r = 1 is often related to the 
radiation brightness temperature. Extending these results to 
Na I Di and Ca II K filtergrams, we are able to use equation]^ 
to interpret out findings with T acting as the radiation bright¬ 
ness temperature. An additional justification for using equa¬ 
tion]^ comes from the fact that due to both filters having non- 
infinitely narrow spectral windows, they will naturally sample 
a considerable portion of the optical continuum. Eor example, 
the Na I Di and Ca II K filters cover 0.17 A and 1 A wide re¬ 
gions, respectively, with continua being clearly visible in the 
form of granulation in the image sequences. 

The inferred phase speeds, which are detailed in Table 1, 
also provide further verification that the waves are slow 
modes. Slow sausage modes are sh own to propagate below 
the sound speed, cs, of the flux tube ( jEdwin & Roberts] 1983[). 
The sound speed inside the pore is derived from equation 
using the atmospheric parameters displayed in Eigure ]2 and 
was estimated as ^7 kms“^ at the photospheric layer, which 
is greater than the estimated phase speeds of the waves in 
each band pass, thus supporting a slow mode interpretation. 
The sound speed estimated in this way is also consistent with 
the numer i cally ca lculated value of the sound speed found in 
Jess et al. (|2012b 1 . Importantly, however, is the fact that the 


phase speeds determined for the atmospheric heights corre¬ 
sponding to the 4170 A continuum and G-band are consistent 
with the phase speeds es tima ted from the time-lag measure¬ 
ments detailed in Section ]3.1] This self-consistency provides 
reassurance in our interpretation that the observed waves are 
upwardly propagating slow axisymmetric oscillations. 

3.4.3. Defining whether the Modes are Surface or Body 

Sausage mode waves can be further categorised based on 
their stracturing within the flux tube. Body modes exhibit pe¬ 
riodic structure along the tube radius for ra dial modes greater 
than z ero, whereas surface modes do not ([Edwin & Roberts] 


1983 I. The distinction between surface and body modes has 
been shown to have an effect on the phase speed of the waves 


where both of the derivatives are to be evaluated at po and 


( Moreels et al.|]20T3 l, but of greater importance are the im- 
plications related to the calculated energy values. As detailed 
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Grant et al. 


Energy Flux of Individual Sausage Mode Waves 
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Height from Solar Surface (km) 

Fig. 5.— The calculated energy flux of 210 s (red) and 290 s (blue) sausage 
mode oscillations, both in W m~^. The hei ghts from the solai* surface are ap¬ 
proximate values that are stated in section [3T| with the energy flux values 
plotted on a logarithmic scale to better highlight the rapid decrease in ob¬ 
served energies with increasing atmospheric height. 


in [Moreels et al. ( |2015b| l, the additional internal oscillations 
associated with hody modes demand that the two cases he 
considered separately. As such, individual energy equations 


for both surface and body modes are defined by Moreels et 
al. (2015b|). Therefore, it is necessary to verify whether the 


oscillations observed in the pore are surface or body modes 
before energy calculations can be undertaken. 

A neat verification of the surface/body nature of sausage 


modes was initially presented by Moreels et al. 
whereby the ratio of wave constants, k, at a given height may 
be defined as. 


mt may 


2 ^ 

icl + vlW4-uj^) ’ 

where P represents the period of an oscillation, k = jPv-p'h 
is the longitudinal wavenumber involving the phase speed of 
the wave, Uph, and ct is the tube speed of wave defined as. 


Ct 


CSVA 

{cl+v\Y 


(7) 


wi th all other terms defi ned in previous sections. It is stated 
in [Moreels et al.| ( |2013| ) that in the case of surface modes, 
the evaluated expression for will be positive. All values 
needed to calculate can be extracted from the model pa¬ 
rameters displayed in Figure or from those defined in Ta¬ 
ble Using these terms, Av was evaluated to be positive 
over all atmospheric heights, indicating that the sausage os¬ 
cillations observed in this pore are slow propagating surface 
modes. 


3.4.4. Calculating the Energy Flux 

With the parameters of the plasma inside and outside the 
pore defined, alongside the phase speeds for individual peri¬ 
odicities at each height and the identification of their surface 
mode nature, the energy flux can sub sequently be calculated. 
The work of [Moreels et al.j ( [2015b[ l utilises a uniform flux 
tube model, and calculates the wave energy by integrating 


over the entire flux cylinder and its exterior, before averag¬ 
ing over the period and wavelength of the oscillation, with 
the equations used in this case being derived for surface mode 
waves. Therefore, the resulting energies depend on the phys¬ 
ical features of the pore (size, magnetic field strength, Alfven 
and sound speeds, and density), in addition to the properties 
intrinsic to the wave (period, phase speed and the Lagrangian 
displacement at the flux tube boundary). For further details, 
a full and detailed description of the theory r equired to cal- 
culate the energy of the waves is presented by [Moreels et al.| 
( [2015bi l. 

The energy flux values for each wave are detailed in Ta¬ 
ble 1, with the energies plotted in Figurej^to highlight the ob¬ 
served dissipation with atmospheric hei^t. It is seen that the 
waves initially contain in excess of 35 000 W m~^ at the so¬ 
lar surface. The estimated wave energy decreases with height, 
with their associated energies at the Ca II K level being three 
orders-of-magnitude less than at their formation at the base of 
the magnetic pore. 

These results conhrm the suitability of a magnetic pore to 
carry compressive wave energy into the chromosphere. The 
observed energy damping with height is directly related to 
the reduction in observed wave power at higher atmospheric 
heights. This may have physical meaning since the structuring 
of the pore is more noticeably diffuse at upper-chromospheric 
heights, hence suggesting the pore becomes unsuitable for ef- 
hcient wave energy transportation. It is also expected that 
wave power will be lost due to the reflection of a portion of 
the waves. Our phase angle histograms suggest that approx¬ 
imately 30% of identihed wave activity between individual 
bandpasses propagates downwards. This can lead to a reduc¬ 
tion in the observed wave energetics at higher atmospheric 
heights and, in part, explain the drop in energy flux from the 
photosphere to the chromosphere. 

An explanation of the sharp drop in energy flux with atmo¬ 
spheric height could be attributed to leaky wave characteris¬ 
tics, i.e., waves that are f reely able to release their energy into 
the surrounding plasma ( Cally|19^ l. Sausage mode oscilla- 
tions displaying a le aky nature have already been observed by 
[Morton et al.[ ( [2012| ), and our suggestion for this is based upon 
the relatively long wavelengths of the observed oscillations. 
The wavelength of t he oscillations can be deduce d using a 
method pr esented by Molowny-Horas et al. ( 1997| l, and later 
applied by [Zaqarashvili et al.|([2007|l, where the wavelength 
can be dehned as, 

A = ^ , (8) 

dip 


where 5s is the distance traversed by the waves and 5p is the 
phase difference of the wave observed between the two dis¬ 
tinct heights. However, when employing this equation one 
must consider a potential 27r uncertainty in the phase angle 
calculated from the phase difference analysis. Such uncertain¬ 
ties are mitigated here as a result of high cadence observations 
that allow us to isolate the start/end points of the oscillation 
cycles in each of the respective bandpasses. As a result, we 
have conhdence that wave activity is observed hrst in hlter- 
grams corresponding to the lower portions of the solar atmo¬ 
sphere, hence validating our interpretation that the waves are 


corresponding to the 4170 A continuum to G-band transition, 
where 5s ^ 75 km and 5p = 6.12°, the wavelength of the 
sausage modes can be estimated as A « 4400 km. This far 
exceeds the distance over which these waves are observed to 


upwardly propagating. Using the values stated in Section 3.1 
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propagate (~ 1000 km). Therefore, it appears that the ob¬ 
served drop in energy flux occurs within one quarter of the 
wavelength of the waves. Many currently proposed damping 
mechanisms hinge on the assumption that multiple cycles of 
the wave are present. As a result, none of these would ade¬ 
quately describe the wave damping we observe in the mag¬ 
netic pore. 

An alternative scenario for the damping of the observed en¬ 
ergy may be through mode conversion. As the waves prop¬ 
agate upwards, they will pass through a region where the 
plasma /? (i.e., the ratio of the plasma pressure to the magnetic 
pressure) equals 1. In this regime favourable conditions ex- 


ist for waves to couple and mode convert (e.g., Ulmschneider 
et al.||199T i. For magnetic regions of the Sun s atmosphere. 


the 13 = 1 layer is expected to be beneath the typ ical heights 
sampled by the Ca II K Alter Pess et ar]|2012a I, thus mode 
conversion of slow sausage mode s to fast modes wit hin in¬ 
clined magnetic fields may occur (jBogdan et al.|2003|l. Such 
mode conversion would likely occur at the inclined edges of 
the pore perimeter, and coupled with the propensity for fast 
modes to be leaky in nature, would help facilitate the rapid 
dissipation of wave energy into the surrounding plasma. 

We must draw attention to the previously described model 
limitations, alongside other features that contribute uncertain- 


ties to the calculated energies. As mentioned in Section 3.4.2 


the model parameters used to describe the pore contain in¬ 
accuracies at different heights in the atmosphere, which will 
have an effect on the derived energy flux values. As men¬ 
tioned previously, the magnetic held strengths used may be 
overestimated at photospheric heights, which will cause a sub¬ 
sequent overestimation of the wave energy flux in the 4170 A 
continuum and G-band data sets. Alongside this, the potential 
underestimation of the pore density at chromospheric heights 
(i.e., through potentially overinflated expansion factors) may 
lead to the energy content of the waves in the Na I Di and 
Ca II K bandpasses being underestimated. Furthermore, the 
uniform flux tube model employed here has several limita¬ 
tions, i.e., it does not consider several physical effects such as 
density stratification, flux tube expansion and/or dissipation. 
However, while the uniform flux tube model does not con¬ 
sider these effects, it is still valid when used as a first order 
approximation, especially considering the pore studied here 
exhibits minimal expansion with atmospheric height. Thus, 
flux tube expansion will not have a significant impact on the 
overall energy flu x since the 15% area exp ansion observed is 
relatively small. |Andries & Cally| (|201 1)1 have shown that, 
for a slowly expanding flux tube, the perpendicular eigen- 
functions (which are used here to calculate the energy flux; 
[Moreels et al.|2015b| l of MHD waves remain almost constant. 
In addition, density stratification is not accounted for in the 
equilibrium model to calculate the energy flux. On the other 
hand, the energy flux is calculated at each localised height 
separately using the equilibrium density at that height. Thus, 
density stratification is to some degree naturally included in 
the energy flux calculations. 

For future studies, combining a more realistic pore model 
atmosphere alongside a more detailed flux tube model may 
subsequently result in improved energy and dissipation es¬ 
timates, particularly for atmospheric heights corresponding 
to the Na I Di and Ca II K image sequences. Furthermore, 
for display purposes in Figure we utilised the Na I Di 
and Call K f ormation heights e s timate d by |Simoniello et 


al. (2008|l and Beebe & Johnson ( 1969|l, which correspond 


to values approximately equal to 400 km and 800 km, respec¬ 
tively. As discussed above, these formation heights are rela¬ 
tively uncertain due to the passbands used sampling portions 
of the spectral line wings, in addition to changes in opacity 
across the magnetically structured pore itself. However, it is 
still observationally evident from Figure[T]that each bandpass 
samples a different region of the solar atmosphere, and as a 
result, still highlights the clear reduction in wave energy with 
atmospheric height even if the specific heights have a degree 
of uncertainty. 


3.5. Oscillations in the Magnetic Field 

In this paper, oscillations in the pore area, intensity and line- 
of-sight velocity have been utilised to investigate the char¬ 


acteristics of s ausage mode waves. However, Fujimura & 
|Tsuneta| ( 2009|l revealed that simultaneous perturbations in 


the associated magnetic flux should also be evident. The au¬ 
thors undertook analyses of the spectropolarimetric proper¬ 
ties of the photospheric Fe I 6302.5 A absorption line, from 
data obtained by the Hinode spacecraft, to detect root mean 
squared oscillation amplitudes between 0.3% and 1.2% above 
the mean. Therefore, an additional investigation of fluctua¬ 
tions in the magnetic held, alongside those already consid¬ 
ered, would provide further information about the mode and 
behaviour of the sausage waves, including the ability to diag¬ 
nose non-trivial situations such as the superposit ion of both 
fast and slow modes wi thin the same structure ( [Moreels 
Van Doorsselaerej 2013|l. Unfortunately, the data presented 


here does not have complimentary high-resolution spectropo¬ 
larimetric observations. However, a natural question arises as 
to whether it is possible to detect magnetic field oscillations, 
particularly for larger-scale structures including solar pores, 
with more coarse resolution instruments such as HMI. Thus, 
the first important step is to calculate what magnetic held am¬ 
plitude would be expected from the sausage mode oscillations 
detected in our present data set, and then compare these to the 
resolution and sensitivity of the HMI instrument. 

The magnetic flux, $, passing through a surface is given by. 


= B • S 


(9) 


where B and S are the vector magnetic held strength and pore 
area, respectively. The vector area is defined as. 


y^^njSj 


( 10 ) 


where h and Si are the normal to the pore area and scalar area 
of the ith element of the pore surface, respectively. Now let us 
consider a surface with a normal vector in the direction of the 
magnetic held (i.e., B = BqZ = B^), providing a pore cross- 
section that is trapped in the x~y plane with a cross-sectional 
area. A, so that hS — zA. Perturbing the variables, linearis¬ 
ing, and assuming the net flux remains constant through the 
surface defined by the cross-sectional area, gives, 

$ = Bo-So + Bi.So + Bo-Si, (11) 


or equivalently, 

0 = b,Ao + BoAi , ( 12 ) 

where Ai and are the perturbations in the area and the ver¬ 
tical (i.e., z direction) magnetic held strength, respectively, 
with Aq and Bq denoting the mean unperturbed area and mag¬ 
netic held values. Perhaps more importantly, we can re-write 
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this relationship as. 


^0 


Bo ■ 


(13) 


Hence, to be able to calculate the magnitude of the change 
in magnetic field strength induced as a result of the upwardly 
propagating sausage-mode waves, we need to utilise the area 
amplitudes derived for each observational bandpass that are 
detailed in Table 1. 

Considering a photospheric average vertical magnetic field 
strength of the pore equal to approximately 900 G (see, e.g., 
the lower-left panel of Figure[T]|, and using the cross-sectional 
area and amplitude values detailed in Table 1, we are able 
to compute the magnitude of bz expected. Given the esti¬ 
mated formation heights for each bandpass, defined in section 
mi the scaled [Maltby et ar] ( |19861 l magnetic field strengths 
corresponding to these heights can be selected (solid line in 
the lower-left panel of Figure Thus, we can estimate the 
magnetic field oscillation ampitudes present in each band¬ 


pass to be GOG, 30G, 10G and 5G, for the 4170A contin¬ 
uum, G-band, Na I Di and Ca II K channels, respectively. 
For the HMI magnetogram formation height of ^300 km, and 
using a «2.2% area amplitude similar to that of G-band and 
Na I Di observations, we estimate that a magnetic field os¬ 
cillation equal to 19 G should be present in the HMI magne¬ 
tograms. Such an oscillation with a peak amplitude «19G 
is close to the noise limit inherent to the HMI instrument. 


where previous analyses by Liu & Schuck|(|2012|l and Welsch 
|& Fisher] ( |2012| l have detected uncertainties in the measured 
field strengths to be wlOMxcm”^ and wdOG, respectively. 
Furthermore, the cadence of the HMI vector (720 s) and line- 
of-sight (45 s) magnetograms are not ideally suited to accu¬ 
rately detect waves with periodicities as short as 3 minutes, 
suggesting the reliable detection of magnetic field oscillations 
in pores is unlikely with the HMI instrument. However, ex¬ 
ploiting this relationship in the future will allow comparisons 
to be made between observed changes in magnetic isosur- 
faces and the oscillat ions found within Stokes profiles (e.g., 
[Balthasar et al.| 2000| ). In addition, as there are known diffi¬ 
culties interpreting complex Stokes measurements, observing 
the changes in pore area could provide a novel way of val¬ 
idating the magnetic oscillations derived from simultaneous 
Stokes profiles. 


4. CONCLUSION 

Observations of sausage mode waves propagating from the 
solar photosphere to the base of the transition region are 
presented here for the first time. High resolution, multi¬ 
wavelength datasets of a magnetic pore are investigated using 
wavelet and Fourier analyses to identify oscillations within 
the pore structure. Simultaneous oscillations in both intensity 
emission and the area encompassed by the pore provide ev¬ 
idence of slow sausage mode activity. Through analysis of 
these oscillations in adjoining bandpasses, it is identified that 
sausage modes with periods between 181 — 412 s are pre¬ 
dominantly propagating upwards from the photosphere to the 
chromosphere. The phase angles calculated for these waves 
between adjacent bandpasses indicate a phase speed in ex¬ 
cess of 15 kms“^. However, the line-of-sight velocity in¬ 
formation for this pore shows a strong simultaneous plasma 
upflow, which is believed to dominate these high speeds. As 
a result, the initial velocity calculation can be considered as 
the Doppler-shifted phase speed of the waves. Using theo¬ 


retical considerations, the waves are verified to be slow sur¬ 
face modes, propagating with a phase speed on the order of 
3 kms“^. The observed periods of the waves are similar to 
that of the global solar acoustic oscillations, and provide fur¬ 
ther insight that sausage mode waves may be driven by the 
ubiquitous p-modes. 

For the first time, a multi-layer analyses of the energy 
content of the sausage mode waves is documented. Indi¬ 
vidual sausage modes were identified for further study, with 
their specific oscillation amplitudes in intensity and area de¬ 
duced. It can be seen in Table[^that the oscillation amplitudes 
increased between photospheric and chromospheric heights, 
consistent with what would be expected as the waves propa¬ 
gate into less dense regions, since there will be less external 
restoring forces for the area perturbations. Through the iden¬ 
tification of these amplitudes, the energy flux of the waves 
at each bandpass was calculated. It is seen that the waves 
are formed at the solar surface with a large energy magni¬ 
tude in excess of 35 000 W m“^, which significantly reduces 
as they propagate into the chromosphere. The limitations 
of the energy flux model are also presented, and we suggest 
how improved atmospheric and flux tube models may be able 
to help more accurately constrain the calculated wave ener¬ 
gies as a function of atmospheric height. However, impor¬ 
tantly we show that magnetic pores can act as viable conduits 
for sausage mode waves to transport energy to higher atmo¬ 
spheric heights. The energy flux values calculated show that 
these waves have lost the vast majority of their energy be¬ 
fore they reach the chromospheric layer observed in Ca II K 
filtergrams. This could be due to wave reflection and conven¬ 
tional damping methods, which are difficult to identify obser- 
vationally, or that these waves are inherent leaky modes and 
freely able to radiate energy into the surroundings. However, 
MHD waves previously observed in chromospheric structures 
have exhibited far greater localised energy fl ux, with Alfvenic 
waves identified to contain 4 — 7 kW m~^ ( De Pontieu et al. 
2007 1 , Alfven waves with energy flux of 15 kW m“^ Pess 


et al.||2009| an d sausage modes di splaying approximately 


To W in“^ ( Morton et al.|201^ . The smaller chromo¬ 
spheric energy values detected in our observations, compared 
to previous studies, could be a result of the structuring of the 
pore itself. The magnetic pore under consideration is small 
and dynamic, is therefore more difficult to resolve from the 
background in chromospheric images, and cannot be posi¬ 
tively identified in simultaneous higher atmospheric imaging 
by SDO. The lack of pore structuring at higher heights may be 
an indication of it’s instability, and thus the pore is unable to 
remain a viable energy conduit for compressible MHD waves 
at these higher atmospheric heights. 

Due to any proposed heating occuring below observable 
MHD scales, there is a lack of observational evidence for the 
heating of chromospheric plasma as a result of the damping 
of sausage mode energy seen in the pore. Instead, the results 
presented here help to clarify the validity of magnetic pores 
as energy conduits for sausage mode waves. Further investi¬ 
gation of sausage modes in a variety of pores at high spatial, 
temporal and spectral resolutions may provide a clearer pic¬ 
ture on the energy content and dissipative potential of these 
waves. With the advent of the Interface R egion Imaging Spec¬ 
trograph (IRIS; De Pontieu et al.|20I4| l mission, and the un¬ 
precedented resolution capabilities of the upcoming Daniel 
K. Inouye Solar Telescope (DKIST, formerly the Advanced 


Technology Solar Telescope, ATST; Keil et akl 2003 Rim- 
jmele et al.|2010|l, there will be new and novel ways to probe 
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the importance that magnetic pores and sausage mode waves 
play in the heating of both the chromosphere and corona. 
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